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We report on iC*" production at mid-rapidity in Au+Au and Cu+Cu collisions at -^snn = 62.4 
and 200 GeV collected by the Solenoid Tracker at RHIC (STAR) detector. The /sT*" is reconstructed 
via the hadronic decays /f*" — > K^tx~ and K*^ — > K~tt'^ . Transverse momentum, pT, spectra are 
measured over a range of pT extending from 0.2 GeV/c up to 5 GeV/c. The center of mass energy 
and system size dependence of the rapidity density, dN/dy, and the average transverse momentum, 
(Pt), are presented. The measured N [K*^) / N [K) and N {(p) / N {K*^) ratios favor the dominance of 
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re-scattering of decay daughters of K*" over the hadronic regeneration for the production. In 
the intermediate px region (2.0 < pT < 4.0 GeV/c), the elliptic flow parameter, V2, and the nuclear 
modification factor, Rcp, agree with the expectations from the quark coalescence model of particle 
production. 
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I. INTRODUCTION 

The main motivation for studying heavy-ion collisions 
at high energy is the study of Quantum Chromodynam- 
ics (QCD) in extreme conditions of high temperature and 
high energy density ^ . Ultra- relativistic nucleus-nucleus 
collisions at the Relativistic Heavy-Ion Collider (RHIC) 
create nuclear matter of high energy density over an ex- 
tended volume, allowing QCD predictions to be tested in 
the laboratory. At high temperature and density, QCD 
predicts a phase transition from nuclear matter to a state 
of deconfined quarks and gluons known as the quark 
gluon plasma (QGP). One of the proposed signatures of 
the QGP state is the modification of vector meson pro- 
duction rates and their in-medium properties 0-131 • 

The K* meson is of particular interest due to its very 
short lifetime and its strange valence quark content. This 
makes the K* meson sensitive to the properties of the 
dense matter and strangeness production from an early 
partonic phase 0. Since the lifetime of the K* is 
~ 4 fm/c, less than the lifetime of the system formed 
in heavy ion collisions 7], the K* is expected to de- 
cay, re-scatter, and regenerate all the way to the kinetic 
freeze-out (vanishing elastic collisions). The character- 
istic properties of the resonance may be modified due 
to high density and/or high temperature of the medium 
causing in-medium effects. Various in-medium effects are 
partonic interaction with the surrounding matter, the in- 
terference between different scattering channels, and ef- 
fects of phase space distortion due to re-scattering of the 
decay daughter particles [1, Q . Measurement of the K* 
meson properties such as mass, width, and yields at var- 
ious transverse momenta can provide insight for under- 
standing the dynamics of the medium created in heavy- 
ion collisions. 

Of particular interest in resonance production is un- 
derstanding the role of re-scattering and regeneration ef- 
fects. Due to the short K* lifetime, the pions and kaons 
from the K* that decay at the chemical freeze-out re- 
scatter with other hadrons. This would then inhibit the 
reconstruction of the parent K*. However, in the pres- 
ence of a large population of pions and kaons, these may 
scatter into a K* resonance state and thus contribute to 
the final measured yield [13] ■ The interplay of these two 
competing processes becomes relevant for determining 
the K* yield in the hadronic medium. These processes 
depend on the time interval between chemical (vanish- 
ing inelastic collisions) and kinetic freeze-out, the source 
size, and the interaction cross section of the daughter 
hadrons. Since the tttt interaction cross section [Tl| is 
larger than the nK interaction cross section jT^, the fi- 



nal observable K* yield may decrease compared to the 
primordial yield. A suppression of the yield ratio such as 
N{K*)/N{K) or N{K*)/N{(j)) is expected in heavy ion 
collisions compared to the same inp+p collisions at simi- 
lar collision energies. This suppression can be used to set 
a lower limit on the time difference between the chemical 
and the kinetic freeze-out [3, [isj . The experimental data 
on the system size, beam energy, and centrality depen- 
dence of this suppression can be used to correlate the life- 
time of the fireball with its size. Although the measured 
values of the resonance yield, mean pT and the elliptic 
anisotropy coefficient V2 are all expected to be affected 
by coUisional dissociation processes, semi-hard scatter- 
ing and jet fragmentation, the measurements presented 
in this paper have been discussed within the framework 
of re-scattering and regeneration. 

The nuclear modification factors such as Raa and 
Rcp are of vital importance in differentiating be- 
tween the effect of hadron mass and hadron type (baryon 
or meson) in the particle production. In the intermediate 
PT range (2.0 < px < 4.0 GeV/c), Rcp of A (a baryon) 
and Kg (a meson), as measured by STAR, are different. 
The observed differences are understood as coming from 
differences in particle type (the baryon-meson effect), in 
agreement with the quark coalescence model Be- 
cause the mass of the K* meson is comparable to the 
mass of the A baryon, it is interesting to compare the 
Rcp of K* with those of Kg and A to check whether the 
results confirm to the expectations of the quark coales- 
cence model. Previous measurements of Rqp for K* in 
Au-f Au collisions at -y/sNJV = 200 GeV were not precise 
enough to make such a conclusion ijj. In this paper we 
present a measurement of Rcp of the K* from a higher 
statistics data set collected in the year 2004. 

In the intermediate px range, the elliptic flow param- 
eter, V2, for different hadrons shows a deviation from 
the particle mass orderin g as seen in the low px regime 
(pT < 1.5 GeV/c) [13, [Bl llTI . For identified hadrons, V2 
follows a scaling with the number of constituent quarks, 
n, as expected from the quark coalescence model P.[l^. 
The K* meson is expected to follow the scaling law with 
n = 2. The K* produced via regeneration of kaons and 
pions during hadronization, on the other hand, would 
follow the n = 4 scaling [l^. Previous STAR measure- 
ments with a smaller data sample found n = 3 ± 2 [^ 
and could not conclusively determine the K* production 
mechanism. The additional V2 data presented in this pa- 
per, for Au-|-Au collisions at 200 GeV, may conclusively 
provide information about the K* production mechanism 
in the intermediate px range. 

In previous STAR measurements, K* production was 
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studied using data from Au+Au, p+p, and d+Au colli- 
sions at 200 GeV d ^ and Au+Au collisions at 130 
GeV [5]. The hadronic decay channels used in these 
analyses were K*° — ^ K+Tr~ , K*^ — >• if~7r+, and 
K*^ — >■ Kg + TT^. In this paper we present new data on 
the pt distribution, (pt), and dN/dy of K*'^ in Au+Au 
collisions at ^snn = 62.4 GeV and Cu+Cu collisions at 
./sWn = 62 A GeV and 200 GeV. The data sample for 
Au+Au collisions at 200 GeV is 6.5 times larger than 
previous measurements, allowing us to make more quan- 
titative conclusions from the V2 and Rcp measurements. 
This broad systematic study, with two different colliding 
beam energies and two different colliding species, enables 
us to study the system size and energy dependence of 
various K*^ properties in heavy ion collisions. To reduce 
statistical errors, the samples of K*'^ and K*^ were com- 
bined and are referred to as K*'^ in the present work, 
unless specified otherwise. 

The paper is organized as follows. In Section II we dis- 
cuss the detectors used in this analysis and details of the 
analysis procedure. For more details on the mixed event 
procedure used to extract the K*'^ yields, the systematic 
uncertainty estimation, and the procedure to obtain W2, 
we refer the reader to our earlier publications @ . Our 
results on pT spectra, dN/dy, (pt), particle ratios, V2, 
and i?cp of K*^ are presented in Section III. The results 
are summarized in Section IV. 



II. EXPERIMENT AND DATA ANALYSIS 

The results reported here, represent data taken from 
Au+Au coUisions at y/sWW = 62 A GeV and 200 GeV in 
the year 2004 and Cu+Cu collisions at ^/snn = 62 A GeV 
and 200 GeV in the year 2005, using the STAR detector 
at RHIC j2Q,] . The primary tracking device within STAR, 
the Time Projection Chamber (TPC) [2l|, was used for 
the track reconstruction of the decay daughters of K*'^. 
The TPC provides particle identification and momentum 
information of the charged particles by measuring their 
ionization energy loss, dE/dx [2ll |. 

The data were collected with a minimum bias (MB) 
trigger. In Au+Au collisions the MB trigger requires 
a coincidence between two zero degree calorimeters 
(ZDC) nil. The ZDCs are located 18 m away from the 
nominal collision point (center of TPC), in the beam di- 
rection, at polar angle, 9, less than 2 mrad. For Cu+Cu 
collisions at 62.4 GeV, the minimum bias trigger was a 
combination of the signals from the ZDC and the Beam 
Beam Counter (BBC). The BBC at 3.3 < r/ < 5.0 com- 
pensates for the trigger inefficiency of the ZDC in cen- 
tral events. To ensure uniform acceptance in the pseudo- 
rapidity, 77, range studied, events with primary vertex 
position, Vz, within ± 30 cm from the center of the TPC 
along the beam line were selected. 

Centrality is defined as function of the fractional cross 
section measured as function of the uncorrected charged 
particle multiplicity within the pseudo-rapidity window 



1 77! < 0.5 for all events @,[2^. The most peripheral events 
were not taken into account due to large trigger and ver- 
tex finding inefficiencies. Table U lists all the collision 
systems studied with the Vz cut, centrality range, and 
number of events used in the analysis. 

TABLE I: List of datasets used in the analysis. Cuts on Vz, 
centrality range selected and number of events used are also 
given. 



Collision systems 


Centrality 


\Vz 


cm 


Events 


Au+Au (62.4 GeV) 


0-80% 


< 


30 


7 X 10" 


Cu+Cu (62.4 GeV) 


0-60% 


< 


30 


8 X 10'^ 


Au+Au (200 GeV) 


0-80% 


< 


30 


13 X 10^ 


Cu+Cu (200 GeV) 


0-60% 


< 


30 


19 X 10'' 



Figure [T] shows the typical dE/dx measured by the 
TPC in Au+Au collisions at 200 GeV in year 2004 as 
a function of momentum, p, divided by charge of the 
particle, q. The different solid lines in the Figure [1] rep- 
resent modified Bethe-Bloch predictions for different par- 
ticle species [23 - |26j . The charged pions and kaons can 
be separated in momenta up to about 0.75 GeV/c while 
(anti-)protons can be separated in momenta up to about 
1.1 GeV/c. The particle identification can be quantita- 
tively described by the variable Na, which for pions is 
defined as: 

7V7- f 1 \dE / dxjrneasured /-, \ 

Nan = p log ^ ^ 1 

R < dE/dx 

where dE / dxmeasured is the measured energy loss for a 
track, {dE/dx)Tr is the expected mean energy loss for a 
pion track at a given momentum [2^ [2^ , and R is the 
dE/dx resolution which is around 8.1%. 




Momemtum/Charge p/q (GeV/c) 

FIG. 1: (Color Online) dE/dx for charged particles versus 
momentum divided by charge of the particle as measured in 
STAR TPC for Au+Au collisions at = 200 GeV. The 

curves are the Bethe-Bloch predictions for different particle 
species. 

K*'^ mesons were reconstructed from their hadronic 
decay channels, K*^ K^tt~ and K*^ if~7r+, us- 
ing charged tracks reconstructed with the TPC. Because 
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the K*'^ decays within a very short time, its daughter 
particles seem to originate from the interaction point. 
Charged kaons and pious with a distance of closest ap- 
proach to the primary vertex (DC A) less than 1.5 cm 
were considered for Au-}-Au collisions at 62.4 GeV and 
Cu+Cu coUisions at 62.4 and 200 GeV. In the case of 200 
GeV Au-f Au collisions, the DCA cut was set at 2.0 cm. 
The charged pion and kaon primary tracks thus selected 
were required to have their respective values less than 
2, with at least 15 fit points inside the TPC. This was 
done to ensure good track fitting with good momentum 
and dE/dx resolution. Further, the ratio of the num- 
ber of fit points to the number of maximum possible fit 
points was required to be greater than 0.55 to avoid se- 
lection of split tracks. In Au-|-Au collisions at 62.4 GeV 
and Cu-hCu collisions at 62.4 GeV and 200 GeV all the 
candidate tracks were required to have jryl < 1 while the 
tracks for the 200 GeV Au-|-Au collisions were required 
to have \rj\ < 0.8 to avoid the acceptance drop at high 
r] range. All tracks selected were also required to satisfy 
the condition that their pT were greater than 0.2 GeV/c. 
All the cuts used for the K*^ analysis are summarized in 
Table m 

TABLE IL List of track cuts for charged kaons and charged 
pions used in the K*^ analysis in Au -I- Au and Cu-|-Cu col- 
lisions at 62.4 GeV and 200 GeV. NFitPnts is the number of 
fit points of a track in the TPC. MaxPnts is the number of 
maximum possible points of the track in the TPC. 



Cut 



Values 



Kaon pt (GeV/c) 
Pion PT (GeV/c) 
NFitPnts 
NFitPnts/Max Puts 
Kaon and Pion 77 

DCA 

Pair Rapidity (y) 



(-2.0, +2.0) 
(-2.0, +2.0) 
(0.2, 10.0) 
(0.2, 10.0) 
> 15 
> 0.55 

\v\ < 1-0 
|?7| < 0.8 (Au+Au 200 GeV) 

< 1.5 cm 
< 2.0 cm (Au+Au 200 GeV) 

\y\ < 0.5 



In a typical event, several hundred tracks originate 
from the primary collision vertex. It is impossible to 
distinguish the tracks corresponding to the decay daugh- 
ters of the from other primary tracks. The K*'^ 
was reconstructed by calculating the invariant mass for 
each unlike-sign Kit pair in an event. The resultant dis- 
tribution consists of the true K*'^ signal and contribu- 
tions arising from random combination of unlike-sign Kn 
pairs. The true K*'^ signal constitutes a very small frac- 
tion of the total invariant mass spectrum. The large ran- 
dom combinatorial background must be subtracted from 
the unlike-sign Kt: invariant mass distribution to extract 
the K*^ yield. This random combinatorial background 
distribution is obtained using the mixed-event technique 
[23I [13, [15. In the mixed event technique, the refer- 
ence background distribution was built with uncorrelated 
unlike-sign Kt: pairs from different events. For generat- 



ing the mixed events, the data sample was divided into 
10 bins in event multiplicity and 10 bins in Vz- Unlike- 
sign Kit pairs from events having similar event multiplic- 
ity and Vz were selected for mixing. This was done to 
ensure that the characteristics of the mixed events gen- 
erated were similar to the actual data. The generated 
mixed event sample was properly normalized to subtract 
the background from the same event unlike-sign invariant 
mass spectrum. The normalization factor was calculated 
by taking the ratio between the number of entries in the 
unlike-sign and the mixed event distributions with invari- 
ant mass greater than 1.2 GeV/c^. The Kit pairs are less 
likely to be correlated in this region. 

Figure [2] shows the background-subtracted and pT- 
integratcd unlike-sign Kt: invariant mass spectra corre- 
sponding to minimum bias Au+Au and Cu+Cu collisions 
at y^snn — 62.4 and 200 GeV. The signal-to-background 
ratio as a function of the Kt: pair for Au+Au col- 
lisions at Y^SNN — 200 GeV is shown in Figure |31 The 
values are of similar order for other collision systems. The 
signal-to-background ratio, S/B, is observed to increase 
with decreasing multiplicity of the events and shows an 
increase with increasing . In the unlike-sign spectrum 
we also could have higher and/or lower Kt: mass reso- 
nant states and non-resonant correlations due to parti- 
cle misidentification and effects from elliptic flow in non- 
central collisions. These effects contribute signiflcantly to 
the residual correlations near the signal @ that are not 
present in the mixed-event sample. These residual corre- 
lations are also subtracted from the background using a 
background function described in Section III . 



III. RESULTS 

A. Mktt peak and width 

The invariant mass distributions (typical distributions 
shown in Figure [5]) for various bins, in Au+Au and 
Cu+Cu collisions at ^snn — 62.4 and 200 GeV were 
flt using a function representing a non-rclativistic Breit- 
Wigner (BW) shape plus a linear residual background 
(RBG). The BW and the RBG parts are as given below. 



BW = 



(Mk^ - Mo)2 + ^ 



RBG = a + bMx^ 



(2) 



(3) 



In the above equations, Mq and Fq are the mass and 
width of the K*'^; a and b are the intercept and slope for 
the linear residual background. 

The variations of Mq and Fq with are shown in 
Figure m The error bars shown correspond to statistical 
uncertainties while the bands represent systematic un- 
certainties. In the low px region (< 1 GeV/c), the mea- 
sured widths are consistent with the Particle Data Group 
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FIG. 2: The Kiv pair invariant mass distribution integrated over the K**^ pr for minimum bias Au+Au and Cu+Cu coUisions 
at y^SNN = 62.4 and 200 GeV after the mixed-event background subtraction. The sohd curve is the signal fit to a Breit-Wigner 
function (Eqn. 2) + Linear function (Eqn. 3) while the dashed line is the linear function representing the residual background. 



(PDG) value of 50.3 MeV/c^ while the measured masses 
are within 2a of the PDG value of 896.0 MeV/c^ (H- 
In the higher pT range (> 1 GeV/c), both the mass and 
width of -fC*° are seen to be consistent with the PDG val- 
ues. We do not observe any significant dependence of A'*° 
mass and width on beam energy and colliding ion species 
studied. The systematic uncertainties on the K*'^ mass 
and width measurement were evaluated bin-by-bin, as a 
function of px^ (1) an uncertainty on the signal fit was 
evaluated by replacing the non-relativistic BW function 
with a relativistic BW function, (2) an uncertainty on 
the background was evaluated by varying residual back- 
ground functions by using higher order polynomials, and 
(3) an uncertainty on track selection was calculated by 
varying the particle identification criteria and different 
cuts on the daughter tracks. In the above analysis, low 
Pt kaon tracks were corrected for energy loss due to mul- 
tiple scattering in the detector [2^ |26| | 



B. Transverse momentum spectra 

The K*'^ invariant yields as a function of pt were eval- 
uated by correcting the extracted raw yields for detec- 
tor acceptance and reconstruction efficiency. The raw 
yield was obtained by fitting the data to the BW-f RBG 
function. The efficiency x acceptance was obtained by 
embedding Monte Carlo (MC) simulations of kaons and 
pions from K*^ decays into the real data using STAR 
GEANT and passing these embedding data through the 
same reconstruction chain as for the real data [1^ . In ad- 
dition, the yields were corrected for collision vertex find- 
ing efficiency and the decay branching ratio of 0.66. The 
vertex finding efficiency is 94.5% for Au-|-Au collisions at 
62.4 GeV and 92.2% for Cu-f Cu colhsions at 62.4 and 200 
GeV. The vertex finding efficiency for Au-|-Au collisions 
at 200 GeV is 100%. The variation of efficiency x ac- 
ceptance with pt, for various centralities in the Au-|-Au 
and Cu-I-Cu system for ^snn = 200 GeV, is depicted in 
Figure [5] The absence of centrality dependence in the 
efficiency x acceptance for the Cu-|-Cu system is due to 
small variation in total multiplicity across the collision 
centrality studied compared to those for the Au-f An sys- 
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FIG. 3: The signal-to-background ratio for K* measure- 
ments as a function of pt for different collision centrality bins 
in Au-I-Au collisions at 200 GeV. 
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tern. 

Figure ini shows the pt spectra of K*'^ at mid-rapidity 
(|y| < 0.5) in Au-|-Au and Cu-|-Cu coUisions at ^snn = 
62.4 and 200 GeV for different collision centralities. The 
dashed lines are the exponential fits to the data. 
The fitting exponential function is defined as 



1 



dN/dy 



^-(mT-Mo)/T 



2T:mT dydrriT 2tiT{Mq + T) 



(4) 



FIG. 4: K**^ (a) mass and (b) width as a function of pt 
for minimum bias Au+Au and Cu+Cu collisions at ^/snn 
= 62.4 and 200 GeV . The dashed line represents the PDG 
values of 896.0 MeV/c^ and 50.3 MeV/c^ for mass and width 
respectively. 



where the inverse slope parameter T and yield dN/dy are 
free parameters. Mq is the mass of the K*^. The above 
function is found to provide good fits to the data for both 
coUision systems. The (pt), obtained using the above 
functional form for the distributions, are presented 
in the following section together with the mid-rapidity 
yields dN/dy. 
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FIG. 6: Mid-rapidity K*'^ pt spectra for various collision centrality bins in Au+Au and Cu+Cu collisions at -^^snn = 62.4 
and 200 GeV. The dashed lines represent the exponential fit to data. The errors shown are quadratic sums of statistical and 
systematic uncertainties. 



C. dN/dy and (pt) 

The dN/dy yield at mid-rapidity plotted as a func- 
tion of average number of participating nucleons, (iVpart) , 
is shown in Figure [71 The dN/ dy for K*'^ presented here 
was calculated by using the data points in the measured 
range of the px spectrum while assuming an exponential 



behavior outside the fiducial range. The K*° integrated 
yield is higher for center of mass energies of 200 GeV than 
62.4 GeV. For collisions at a given beam energy with sim- 
ilar (A'part), the dN/dy is similar for Au-|-Au and Cu-f Cu 
systems. A similar behavior was observed for (j) mesons 
at RHIC [13. 

The ivr*° (pt) at mid-rapidity plotted as a function of 
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TABLE III: The K*° dN/dy and (pt) at \y\ < 0.5 measured in Au+Au and Cu+Cu collisions at 62.4 GeV and 200 GeV for 
different collision centralities. Both statistical and systematic uncertainties are given. The error given for T is statistical only. 



Collision systems Centrality (Npgrt) dN/dy T (GeV) (pr) (GeV) 



Au+Au(62.4 GeV) 


0-20% 
20-40% 
40-60% 
60-80% 


275 
137 
60 
19 


6.4 ± 0.4 ± 0.7 
2.9 ±0.2 ±0.4 
1.4 ±0.1 ±0.2 
0.56 ±0.03 ±0.07 


0.36 ±0.02 
0.35 ±0.02 
0.32 ±0.01 
0.26 ±0.02 


0.93 ±0.05 ±0.09 
0.92 ± 0.05 ± 0.08 
0.86 ± 0.06 ± 0.08 
0.75 ± 0.04 ± 0.07 


Cu+Cu(62.4 GeV) 


0-20% 
20-40% 
40-60% 


84 
44 
20 


2.07 ±0.07 ±0.30 
1.15 ±0.06 ±0.20 
0.51 ±0.03 ±0.07 


0.35 ±0.01 
0.33 ±0.01 
0.31 ±0.01 


0.94 ±0.05 ±0.06 
0.86 ± 0.04 ± 0.06 
0.83 ± 0.04 ± 0.06 


Au+Au(200 GeV) 


0-10% 
10-40% 
40-60% 
60-80% 


326 
173 
62 
20 


9.05 ±0.57 ±1.01 
5.43 ±0.17 ±0.69 
1.95 ±0.07 ±0.18 
0.53 ±0.03 ±0.08 


0.41 ±0.02 
0.43 ± 0.02 
0.39 ±0.02 
0.40 ± 0.03 


1.09 ±0.06 ±0.094 
1.09 ±0.03 ±0.12 
1.02 ±0.03 ±0.09 
1.05 ±0.05 ±0.13 


Cu+Cu(200 GeV) 


0-20% 
20-40% 
40-60% 


86 
46 
21 


2.96 ±0.12 ±0.30 
1.55 ±0.06 ±0.20 
0.73 ± 0.03 ± 0.09 


0.40 ±0.01 
0.40 ± 0.02 
0.40 ± 0.02 


1.0 ±0.04 ±0.08 
1.0 ±0.04 ±0.09 
0.95 ±0.02 ±0.07 


p+p (200 GeV) 




2 


0.005 ± 0.002 ± 0.006 


0.20 ±0.01 


0.81 ±0.02 ±0.14 



TABLE IV: The contributions for various sources for estimating the total systematic uncertainties for iT*" at mid-rapidity 
{\y\ <0.5) on dN/dy and (pr) in 0-20% Au±Au collisions at 62.4 GeV. The systematic uncertainties are similar for other 
collision systems. 



Different Sources 


dN/dy (pt) (GeV/c) 


exponential fit 


0% 


0% 


Levy fit 


5.24% 


0.44% 


Background function 


4.49% 


2.88% 


(higher order polynomial) 






Relativistic Breit Wigner 


2.64% 


0.542% 


\Vz\ < 20 cm 


1.6% 


1.05% 


Track type (/C*°) 


3.6% 


4.3% 


Track type (^*° ) 


4.05% 


5.96% 


NFitPnts = 22 


4.32% 


1.35% 


\N^n,N^K\ < 3 


4.45% 


4.87% 


Total Sys. Uncertainty 


11.18% 


9.47% 



This indicates that the (px) is strongly coupled with the 
mass of the particle, in agreement with similar observa- 
tions made previously in Au+Au and d+Au collisions 
at 200 GeV gUi. Table Im] Hsts the dN/dy and (px) 
of K*'^ for various collision systems at different collision 
centralities and beam energies studied. 



The systematic uncertainties on K*'^ dN/dy and (px) 
were estimated as follows j32j : (1) an uncertainty on 
the K*'^ signal fit of the invariant mass spectrum was 
evaluated by replacing the non-relativistic BW function 
with a relativistic BW function, (2) an uncertainty on 
the background distribution fit was evaluated by using 
a higher order polynomial function, (3) by varying the 
track types {K*'^ and K*^\ (4) using different functions 
such as a Levy function |33l , |34| to fit the spectra, (5) an 
uncertainty on track selection was estimated by varying 
the track cuts such as N^ cut, NFitPnts cut, and (6) by 
varying the Vz cut from 30 cm to 20 cm. The systematic 
uncertainties coming from the different sources are listed 
in Table |lV]for Au±Au colhsions at 62.4 GeV. 



(-/Vpart), is shown in Figure[8]for Au+Au and Cu+Cu col- 
lisions at ^/iiS" = 62.4 and 200 GeV. The (px) for K*° 
presented here was calculated by using the data points in 
the measured range of the px spectrum while assuming 
an exponential behavior outside the fiducial range. No 
significant centrality and colliding ion size dependence 
could be observed. However, the (px) values for collisions 
at 200 GeV are seen to be slightly higher than those from 
62.4 GeV in both Au+Au and Cu+Cu collisions. Previ- 
ous measurements of (px) of K*^ in heavy ion collisions 
have been shown to be higher than the corresponding val- 
ues in p+p collisions [6*1. This may be understood from 
the following. According to UrQMD transport model cal- 
culations [30| , K*'^s are more likely to be reconstructable 
in the high px region than in the low px region. This 
is because high px K*^s are more likely to escape the 
medium before the kinetic freeze-out stag e and arc thus 
less affected by in-medium effects [Sll, |32| 

Figure |9] shows the (px) of different particle species (tt, 
K, p, and K*°) in Au+Au collision at ^snn = 62.4 GeV 
as a function of (iVpart)- The (px) of K*'^ is higher than 
the (px) of kaons and pions and closer to that of protons. 
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of the average number of participating nucleons, (A^part), for 
Au+Au and Cu+Cu coUisions at y/swN = 62.4 and 200 GeV. 
The boxes represent the systematic uncertainties. 
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FIG. 8: The mid-rapidity K'^ {pt) as a function {Np^n) for 
Au+Au and Cu+Cu collisions at -^snn = 62.4 and 200 GeV. 
The boxes represent the systematic uncertainties. 



D. Particle Ratio 

Figure [TOT a) shows the ratio of K*'^ and K~ yields, 
N{K*")/N{K-) as a function of (TVpart) in Au+Au and 
Cu+Cu collisions at y/suN = 62.4 and 200 GeV. From 
the figure, no clear beam energy or system size depen- 
dence is observed. The FigfTOlb) shows the K*'^/K~ ra- 
tio in Au+Au, Cu+Cu, and d+Au collisions normalized 
by their corresponding values measured in p+p collisions 
at yiiiF = 200 GeV. This N{K*°)/N{K~) double ratio 
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The mid-rapidity (pt) of tt, K, p and K*^ as a func- 
(Npait) for Au+Au collisions at -^snn = 62.4 GeV. 



is seen to be much smaller than unity in central Au+Au 
collisions. In contrast the N{K*°)/N{K~) double ra- 
tio is close to unity for d+Au collisions. This suggests 
strong re-scattering of decay daughters of K*^ meson, 
resulting in the loss of reconstructed K*^ signal. The 
re-scattering of if*" daughter particles depends on cr^^ 
which is considerably larger than Uj^k, but a.,:K is re- 
sponsible for regeneration of the K*^ meson. Therefore, 
we expect a decrease of the N{K*^)/N{K~) yield ratio 
in heavy ion collisions owing to strong re-scattering of 
K*^ daughter particles. The observed decrease in the 
K*^ /K" double ratio indicates an extended lifetime for 
the hadronic phase as we move from p+p and d+Au to 
Au+Au collisions. The extended lifetime enhances the 
re-scattering effect. Figure [TOTc) shows the energy de- 
pendence of the N{K*^)/N{K^ ) ratio for minimum bias 
Au+Au and Cu+Cu collisions at ^snn = 62.4 and 200 
GeV. Also included in the figure are values obtained from 
p+p collisions at 63 GeV [sl and 200 GeV fe]. At both 
energies, the N{K*'^)/N{K~) for p+p collisions is higher 
than the values in the heavy ion collisions. This can be 
attributed to larger re-scattering of K*'^ daughter parti- 
cles in heavy ion collisions. 

Another ratio of considerable interest is the 
N{(j))/N{K*°) ratio as both the (j) and K*° have the 
same spin and similar mass, but different strangeness 
and lifetime. The lifetime of the cj) meson is 40 fm/c 
(^10 times that of K*^). Due to the relativel y lo nger 
lifetime of the (j) meson and negligible aKK [23|, we 
expect both the re-scattering and regeneration effects to 
be negligible for the (j) meson. Since 4> has two strange 
quarks and K*'^ has one, N {(j)) / N {K*^) can also give 
information regarding strangeness enhancement. 

Figureinia) depicts the N {(j)) / N {K*°) ratio as a func- 
tion of (A^part), corresponding to Au+Au and Cu+Cu 
collisions at ^/snn = 62.4 and 200 GeV. We observe that 
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TABLE V: The mid-rapidity N{K*'^)/N{K') and N{4>)/N{K*°) yield ratio in Au+Au and Cu+Cu collisions at = 62.4 

GeV and 200 GeV for different centralities. The first uncertainty is statistical and the second one is systematic. 

Collision systems Centrality jV(A'*")/jV(A'- ) N{(j>)/N{K*") 
Au+Au(62.4 GeV) 0-20% 0.24 ± 0.02 ± 0.04 0.51 ± 0.03 ± 0.09 
20-40% 0.26 ± 0.02 ± 0.04 0.49 ± 0.03 ± 0.08 
40-60% 0.30 ± 0.03 ± 0.04 0.39 ± 0.03 ± 0.06 

60-80% 0.44 ± 0.04 ± 0.07 0.24 ± 0.02 ± 0.04 

Cu-fCu(62.4 GeV) 0^20% 0.29 ± 0.01 ± 0.05 0.51 ± 0.02 ± 0.1 
20-40% 0.34 ± 0.02 ± 0.04 0.43 ± 0.02 ± 0.08 

40-60% 0.36 ± 0.02 ± 0.05 0.37 ± 0.02 ± 0.07 

Au-FAu(200 GeV) 0-10% 0.20 ± 0.03 ± 0.03 0.82 ± 0.05 ± 0.08 
10-40% 0.26 ± 0.02 ± 0.04 0.68 ± 0.02 ± 0.08 
40-60% 0.29 ± 0.03 ± 0.04 0.58 ± 0.02 ± 0.05 

60-80% 0.25 ±0.03 ±0.04 0.62 ± 0.05 ± 0.1 

Cu±Cu(200 GeV) 0-20% 0.27 ± 0.01 ± 0.03 0.66 ± 0.03 ± 0.09 
20-40% 0.30 ± 0.01 ± 0.03 0.58 ± 0.02 ± 0.09 

40-60% 0.33 ± 0.01 ± 0.04 0.45 ± 0.02 ± 0.09 

p±p(200 GeV) MB 0.34 ± 0.01 ± 0.05 0.36 ± 0.02 ± 0.07 



the ratio tends to increase with increasing (A^part) at a- 
given beam energy. The N {(j)) / N [K*^) ratio is higher 
for y'SNN = 200 GeV compared to ^snn = 62.4 GeV 
for the various collision centralities. At a given beam 
energy and (A'^part) the N{(l)) IN{K*°) ratio is similar for 
Au+Au and Cu+Cu collisions. Figure fllTb) also shows 
N{(f>)/N{K*") ratio in Au+Au, Cu+Cu, and d+Au colli- 
sions normalized by their corresponding values measured 
in p+p collisions at ^snn = 200 GeV. We observe that 
this double ratio increases with collision centrality, which 
favors the re-scattering scenario of K*^ daughter par- 
ticles. The observed increase can also have contribu- 
tions from strangeness enhancement in more central col- 
lisions [1^. Figure [TlTc) shows the energy dependence of 
the N {(j>) / N {K*^) ratio for minimum bias Au+Au and 
Cu+Cu collisions at ^snn = 62.4 and 200 GeV and its 
values from p+p collisions. The -^/s = 63 GeV value for 
p+p collisions is from ISR measurements [s^. At both 
the energies the N {</))/ N {K*'^) ratios for p+p collisions 
are lower than the corresponding values in Au+Au colli- 
sions. Furthermore, there is an indication of an increase 
of the value with beam energy. The study of both ra- 
tios N{K*°)/N{K-) and N{(j))/N{K*°), as a function 
of colliding species, collision centrality, and beam energy 
favors the re-scattering scenario over regeneration. 
The values of the ratios along with the associated uncer- 
tainties are shown in Table [V] 



in the same event to generate an event plane subtracted 
azimuthal distribution in the variable $ — {(p ~ ^J2). 
The corresponding distribution, (PN/dprd^, in the az- 
imuthal angle ^ for all the events in a given pt bin were 
then fitted with a function ^[1 + 2tif'"' cos(2<I))] where A 
is a constant. The fitted value of wf'"' was then divided 
by the reaction plane resolution factor to obtain V2 for 
the Pt window considered [37| . 

Figure shows V2 of K*^ as a function of pt in mini- 
mum bias Au+Au collisions at ^snn = 200 GeV. We fit 
the data using the function 

, , an 
V2(PT,n)^— — TTTT-dn (5 

1 + exp[-{pT/n - b)/c) 

where a, b, c, and d are the parameters extracted from 
such a fit to data obtained earlier for Kg and A . 
Here n, the number of constituent quarks, is the only free 
parameter. The best fit of the K*'^ data with the function 
as given in EqnlS] yields a value of n = 2.0 ± 0.4 {x^ /ndf 
= 2/6). A similar fit of the combined results of Run II 
and Run IV data, taken in the years 2002 and 2004, re- 
spectively, also yields an identical value of n = 2.0 ± 0.4 
{x^/ndf = 4/10). This indicates that K*° are domi- 
nantly produced from direct quark combinations, and the 
regenerated K*'^ component in the hadronic stage is neg- 
ligible compared to the primordial K*'^. 



E. Elliptic Flow 



We apply the standard reaction plane method as em- 
ployed in Refs [s^, [s^l for the analysis of elliptic flow. 
Here, for a given px window, the second order reaction 
plane angle, V'2, was determined event- by- event. We have 
different event planes for every candidate, because 
for every if*" candidate, its daughter particles are ex- 
cluded from the event plane determination. This was 
later subtracted from the azimuthal angle of each track 



F. Nuclear Modification Factor 

Through a measurement of the nuclear modification 
factors Rcp and Raa, one probes the dynamics of par- 
ticle production during hadronization and in-medium ef- 
fects [H, nil- The nuclear modification factor Rcp, 
which is the ratio of the invariant yields for central to 
peripheral collisions, normalized by number of binary col- 
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FIG. 10: (a) Mid-rapidity N(it"*'')/N(A"-) ratio for Au+Au 
and Cu+Cu collisions at ^snn = 62.4 and 200 GeV as a func- 
tion of (iVpart). (b) Mid-rapidity N{K*°) /N {K') in Au+Au, 
Cu+Cu and d+Au collisions divided by N{K*°)/N{K~) ra- 
tio in p+p collisions at y/swN = 200 GeV as a function of 
(iVpart). (c) Mid-rapidity N{K*°)/N{K-) ratio in minimum 
bias Au+Au, Cu+Cu, p+p collisions as a function of -^snn- 
The boxes represents systematic uncertainties. The value of 
N{K''°)/N{K-) ratio in p+p at 63 GeV is from ISR 
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FIG. 11: (a) Mid-rapidity N {<j)) / N {K*'^) ratio for Au+Au 
and Cu+Cu collisions at -^snn = 62.4 and 200 GeV as a func- 
tion of (A^part). (b) Mid-rapidity N {<t)) / N {K'^) in Au+Au, 
Cu+Cu and d+Au collisions divided by N{(j))/N{K*°) ra- 
tio in p+p collisions at ^snn ~ 200 GeV as a function of 
(iVpart). (c) Mid-rapidity N[cj>) /N{K*^^) ratio in minimum 
bias Au+Au, Cu+Cu, p+p collisions as a function of y/swN- 
The boxes represents systematic uncertainties. The values of 
N{4>)IN{K'°) ratio in p+p at 63 GeV is from ISR iSSj. 
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lisions, Nfjin, is defined as 



Rcp 



[dN/[NundpT)] 



central 



[dN/iNbindpT)]'"'"^^'''''''- 



(6) 



where A^bm is calculated from the Glauber model 26]. 
We expect Rcp to be unity at high px (> 2 GeV/c) 
if nucleus-nucleus collisions were mere superpositions of 
nucleon-nucleon collisions. Any deviation observed from 
unity would indicate the presence of in-medium effects. 
Above = 2 GeV/c, the Rcp of tt^, p + p, Kg, and 
A, as measured by STAR, are found to be significantly 
lower than unity. This suggests a suppression of par- 
ticle production at hig h in central collisions relative 
to peripheral ones [l4 , [l^ 24 1 • Theoretically, this is at- 
tributed to the energy loss of highly energetic partons 
while traversing through the dense medium created in 
heavy ion collisions. We also observe that the Rcp of 
Kg and A are different. Since the mass of K*^ is close to 
that of baryons such as p and A, a comparison of Rcp of 
if*" with those for K^ and A can be used to understand 
whether the observed differences in the Rcp of the Kg 
and the A are tied to the particle mass or the baryon- 
meson effect (T^ |. 

Figure [T^ shows the K*'^ Rcp as a function of px com- 
pared to those for A and Kg [l3|. The shaded band 
around the data points represents the systematic uncer- 
tainties and the band around 1 on the right corner rep- 
resents the normalization uncertainty. For Au-{-Au col- 
lisions at 200 GeV the K*'^ Rcp was obtained from the 
Pt spectra of top 10% and 60-80% centrality classes. For 
Au-|-Au collisions at 62.4 GeV the pr spectra of the top 
20% and 60-80% centrality classes were considered. The 
A and Kg Rc p correspond to the pt spectra of the top 
5% and 60-80% Au-f Au coUisions at 200 GeV [ll|. For 
Pt < 1.8 GeV/c, the Rcp of K*^ in Au-|-Au colhsions at 
200 GeV and 62.4 GeV are smaller than that of A and 



K^. This is consistent with the assumption that the re- 
scattering effect dominates over the regeneration effect 
for K*° at low px- For Au-|-Au cohisions at 200 GeV, 
for px > 1-8 GeV/c, the Rcp of K*^ is closer to that for 
K^ (differing from that of A). Since the masses of the A 
and K*^ are similar, the observed difference seems to be 
due to other than mass. The observed differences might 
arise because the A is a baryon whereas if*" is a meson. 
This supports the quark coalescence picture of particle 
production in the intermediate pr range. 



IV. SUMMARY 

STAR has measured the if*" resonance production at 
mid-rapidity in Au+Au and Cu+Cu collision systems 
at y^SNN = 62.4 GeV and 200 GeV. A large sample of 
Au+Au collision data at ^snn — 200 GeV enables us to 
extend the measurements to px ^ 5 GeV/c. The mea- 
sured dN/dy and (px) of if*" are higher at ^snn = 200 
GeV compared to the corresponding values at ^snn — 
62 A GeV. For a given beam energy, the dN/dy and (px) 
are similar for Au+Au and Cu+Cu collisions at a given 
(A^part)- For ySiii" = 62.4 and 200 GeV the if*" (px) is 
comparable to the same for protons, indicating that the 
(px) trends are dependent on the mass. 

The N{K*°)/N{K-) ratio in central Au+Au coUisions 
at both 62.4 and 200 GeV is much smaller compared 
to the respective values in p+p collisions. The data in- 
dicate that heavy ion collisions provide an environment 
with stronger re-scattering of if*" daughter particles rel- 
ative to regeneration. The N{(f>)/N{K*'^) ratio in central 
Au+Au collisions at both 62.4 and 200 GeV is larger than 
that of p+p collisions, again supporting the dominance of 
re-scattering effects. The increase in the N {(f) / N [K*^) 
ratio as a function of beam energy and collision central- 
ity also suggests strangeness enhancement in heavy-ion 
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collisions. 

The large sample of Au+Au collision data at 200 GeV 
allow for a quantitative estimation of elliptic flow of K*^ 
and the interpretation of the V2 in terms of a scaling 
based on the number of constituent quarks. The results 
support the quark coalescence model of particle produc- 
tion. More explicitly, are dominantly produced from 
direct quark combinations, with a negligible regenerated 
component. At low px, the nuclear modification factor 
for is seen to be similar for Au+Au collisions both 
at 62.4 and 200 GeV. At lower px, Rcp for Au+Au col- 
lisions at ^snn = 200 GcV is lower than that for A and 
Kg, which is consistent with the observation that the re- 
scattering effect dominates over regeneration effect. For 
Pt > 1.8 GeV/c, the K*° Rcp in Au+Au collision at 200 
GeV more closely follows that for Kg, and at the same 
time differs from that for A. This also provides support 
for the quark coalescence picture at the intermediate pt 
ranges studied. 
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